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Abstract 

We study the phenomenon of system size stochastic resonance within the nonequilibrium po- 
tential's framework. We analyze three different cases of spatially extended systems, exploiting the 
knowledge of their nonequilibrium potential, showing that through the analysis of that potential 
we can obtain a clear physical interpretation of this phenomenon in wide classes of extended sys- 
tems. Depending on the characteristics of the system, the phenomenon results to be associated 
to a breaking of the symmetry of the nonequilibrium potential or to a deepening of the potential 
minima yielding an effective scaling of the noise intensity with the system size. 

PACS numbers: 05.45.-a, 05.40. Ca, 82.40. Ck 



*Electronic address: izus@mdp.edu. ar 
^Electronic address: wio@ifca.unican.es 



1 



I. INTRODUCTION 



The phenomenon of stochastic resonance (SR) — namely, the enhancement of the output 
signal-to-noise ratio (SNR) caused by injection of an optimal amount of noise into a non- 
linear system — configures a counterintuitive cooperative effect arising from the interplay 
between deterministic and random dynamics in a nonlinear system. The broad range of 
phenomena for which this mechanism can offer an explanation can be appreciated in Ref. Q] 
and references therein, where we can scan the state of the art. 

Most of the phenomena that could possibly be described within a SR framework occur 
in extended systems: for example, diverse experiments were carried out to explore the role 
of SR in sensory and other biological functions |2J or in chemical systems These were, 
together with the possible technological applications, the main motivation to many recent 
studies showing the possibility of achieving an enhancement of the system response by means 
of the coupling of several units in what conforms an extended medium 
analyzing the possibility of making the system response less dependent on a fine tuning of 
the noise intensity, as well as different ways to control the phenomenon Q, ^] . 



In previous papers 



121 ] we have studied the stochastic resonant phenomenon in 



extended systems for the transition between two different patterns, exploiting the concept 
of nonequilibrium potential (NEP) HQ- 

This potential is a special Lyapunov functional 
of the associated deterministic system which for nonequilibrium systems plays a role similar 
to that played by a thermodynamic potential in equilibrium thermodynamics Q|. Such a 
nonequilibrium potential, closely related to the solution of the time independent Fokker- 
Planck equation of the system, characterizes the global properties of the dynamics: that is 
attractors, relative (or nonlinear) stability of these attractors, height of the barriers separat- 
ing attraction basins, and in addition it allows us to evaluate the transition rates among the 



different attractors 



13, 



14j . In another recent paper we have explored the characteristics 



of this SR phenomenon in an extended system composed by an ensemble of noise-induced 
nonlinear oscillators coupled by a nonhomogeneous, density dependent diffusion, externally 
forced and perturbed by a multiplicative noise, that shows an effective noise induced bistable 
dynamics |l5[ . The stochastic resonance between the attractors of the noise- induced dynam- 
ics was theoretically investigated in terms of a two-state approximation. It was shown that 
the knowledge of the exact NEP allowed us to completely analyzed the behavior of the 
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output SNR. 
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17| have shown that 



Recent studies on biological models of the Hodgkin-Huxley type 
ion concentrations along biological cell membranes present intrinsic SR-like phenomena as 
the number of ion channels is varied. A related result shows that even in the absence of 
external forcing, the regularity of the collective firing of a set of coupled excitable FitzHugh- 
Nagumo units results optimal for a given value of the number of elements. From a physical 
system point of view, the same phenomenon -that has been called system size stochastic 
resonance (SSSR)- has also been found in an Ising model as well as in a set of globally 
coupled units described by a 4 theory |lj| . It was even shown to arise in opinion formation 
models Q- 

The SSSR phenomenon occurs in extended systems, hence it is clearly of great interest 
to describe this phenomenon within the NEP framework. More, the NEP offers a general 
framework for the study of the dependence of resonant and other related phenomena on any 
of system's parameters. With such a goal in mind, in a recent paper it was shown that 
SSSR could be analyzed within a NEP framework and that, depending on the system, its 
origin could be essentially traced back to a breaking of the symmetry of such a potential. 
Even those cases discussed in Q| could be described within this same framework, and the 
( "effective" ) scaling of the noise with the system's size could be clearly seen. Here, we discuss 



in more detail the cases analyzed in 



21 1 and present a new interesting one, corresponding to 



the study of SSSR in a system that also shows noise induced patterns, the same one studied 
in Q|. We show that in two of the cases -corresponding to pattern forming systems that 
include only local interactions- the problem could be rewritten in such a way as to present 
a kind of "entrainment" between the symmetry breaking of an "effective" potential together 
with a scaling of the noise intensity with the system size. 

The organization of the paper is as follows. In Section |TI] we focus on a simple reaction- 
diffusion model with a known form of the NEP, that presents SSSR associated to a NEP's 
symmetry breaking. In Section HTT1 we analyze the model of globally coupled nonlinear oscil- 
lators discussed in 19:j, and show that it can also be described within the NEP framework, 
with SSSR arising due to a deepening of the potential wells, or through an "effective" scaling 
of the noise intensity with the system size. We start Section IIVI by briefly reviewing the 
model and the formalism to be used for the case of multiplicative noise. In this case, by 
scaling the NEP with the system size, we show that the system's behavior could be associ- 
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ated to a kind of "entrainment" between the symmetry breaking of an "effective" potential 
and a scaling of the noise intensity with the system size. Finally, we present in Section IS 
some conclusions and perspectives. 



II. A SIMPLE REACTION-DIFFUSION SYSTEM 



Brief Review of the Model 



The specific model we shall focus on in this section, with a known form of the NEP, 



22 



231 ] that, with a piecewise 



22 



corresponds to a one-dimensional, one-component model 

linear form for the reaction term, mimics general bistable reaction-diffusion models 
that is with a cubic like nonlinear reaction term. In particular we will exploit some of the 
results on the influence of general boundary conditionsfcalled albedo) found in [2^] as well 
as previous studies of the NEP 14] and of SR U, U B\. 

nun 

The particular non-dimensional form of the model that we work with is [a, 13, |24| 

d 2 




dt' 



D 



dy 2 



(1) 



We consider here a class of stationary structures 4>{y) in the bounded domain y E [— L, L] 
with albedo boundary conditions at both ends, 



rf0 
dy 



Tk(f>{±L), 



y=±L 



where k > is the albedo parameter. It is worth noting that for fc^Owe recover the usual 



case of Neumann boundary conditions (i.e. -$- =0), while for k — > oo what results is 

y y=±L 

the usual Dirichlet boundary conditions (0(±L) = 0). 

Those stationary structures are the spatially symmetric (stable) solutions to Eq. (0) al- 



ready studied in 



24j . The explicit form of these stationary patterns is (see |24| for details) 



<P(y) 



S wh(y c /VD) p'(k, (L + y)/VD) p(k, L/VD)-\ 
1 - cosh(y/y/D) p(k, (L - y e )/VD) p(k, L/^)~\ 
smh(y c /VD) p'(k, (L - y)/VD) p(k, L/VD)-\ 



if ~ L<y < -y c , 
ti-y c <y < y e , 
ify c <y < L, 



(2) 



with p(k, C) = sinh(^) + k cosh(^), and p'(k, Q 



dp g k ^ . The double- valued coordinate y c , 
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Y (arb. units) 



FIG. 1: Typical form of the patterns (j)(y), for D = 1., L = 2, k = 2 and (fic/tfih = 0.193. The 
continuous line corresponds to the stable nonhomogeneous pattern, the dashed one is the unstable 
pattern. In addition we also have the always present, stable, null pattern. 



at which 



is given by 24] 



L -- L -^ 

2 2 



zp{k, L/y/D) ± y z 2 p(k, L/y/D) 2 + 1 - k 2 
1 + k 



(3) 



with z — 1 — 2<ft c /<f)f l (—1 < z < 1). Typical forms of the patterns are shown in Fig. 1. 

When y^- exists and yf < L, this pair of solutions represents a structure with a central 
"excited" zone (0 > <f> c ) and two lateral "resting" regions (0 < C ). For each parameter set, 
there are two stationary solutions, given by the two values of y c . Figure 5 in [24], that we do 
not reproduce here, depicts the curves corresponding to the relation y c /L vs. k, for various 
values of 4> c /4>h- 



Through a linear stability analysis it has been shown 24j that the structure with the 
smallest "excited" region (with y c = y+, denoted by <f> u (y)) is unstable, whereas the other 
one (with y c = y~ denoted by <fii(y)) is linearly stable. The trivial homogeneous solution 
4>o(y) = (denoted by 0o) exists for any parameter set and is always linearly stable. These 
two linearly stable solutions are the only stable stationary structures under the given albedo 
boundary conditions. We will concentrate on the region of values of z, L and k, where both 
nonhomogeneous structures exist. 
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□ 



For the system with the albedo b.c. that we are considering here, the NEP reads 
Fi<P,k,L] = J {-J h0' + W -0 C )] ^+2" (^(y.*)J Hy+-#y,t) 2 . 

Strictly speaking, this is the system's Lyapunov functional, as we are still considering the 
deterministic case. However, in what follows we will always refer to the NEP both, for the 
deterministic and stochastic cases. This functional fulfills the "potential" condition 

where s ^ y ^ indicates a functional derivative. 

Replacing the explicit forms of the stationary nonhomogeneous solutions (Eq.(J2J)), we 
obtain the explicit expression J(| [3| 



F± = ^[0 Uil , k, L] = -<gyt z + <Pl smh(yt/VD) P -^%- jf ] J^\ (6) 



p(k,L/VD) 

while for the homogeneous trivial solution O — 0, we have instead J-[4>q, k, L] = JF° = 0. 

Figure 2 depicts the nonequilibrium potential !F[(f), k, L] as a function of the system size 
L, for a fixed albedo parameter k, and a fixed value of c /0/i (that is, with fixed value of 
z). The curves correspond to the NEP evaluated on the nonhomogeneous structures, J r± , 
whereas the horizontal line stands for JF°, the NEP of the trivial solution. We have focused 
on the bistable zone, the upper branch being the NEP of the unstable structure, where T 
attains a maximum, while in the lower branch (for <fi = <f) Q or <p = <pi), the NEP has local 
minima. We see that when L becomes small, the difference between the NEP for the states 
4> u (y) and 0i (y) reduces until, for L rs 0.72, they coalesce and, for even lower values of L, 
disappear (inverse saddle-node bifurcation). 

It is important to note that, since the NEP for the unstable solution <p u is always positive 
and, for the stable nonhomogeneous structure 0i, T < for L large enough, and JF > for 
small values of L, the NEP for this structure vanishes for an intermediate value L = L* of 
the system size. At that point, the stable nonhomogeneous structure <fii(y) and the trivial 
solution 0o (y) exchange their relative stability. 

For completeness and latter use, in Fig. 3 we show J-[4>, k, L] but now as a function of k, 
for a fixed value of L and the same value of z. Here we see that the initial large difference 
between the NEP for the states 4> u (y) and 0i(y) reduces for increasing k until, for k — > oo, 
the values for Dirichlet b.c. are asymptotically reached. 



FIG. 2: T\(j),k,L\ vs. L: NEP evaluated at the stationary solutions 4>o(y), 4>\{y) and (p u (y)- Here 
k = 3, D = 1., and c /0 h = 0.193. 




FIG. 3: .F[<^, fc, L] vs. fc: NEP evaluated at the stationary solutions <f>o(y), <f>i(y) and 4> u (y)- Here 
L = 1.2, D = 1., and <J) c /<p h = 0.193. 

B. System Size Stochastic Resonance 

In order to account for the effect of fluctuations, we include in the time-evolution equation 
of o„ r .node, (Eq.JIJ) a fluctuation tenn, that we m odel as an add.t.ve noise 
yielding a stochastic partial differential equation for the random field 4>(y,t) 

-4>{y, t) = D —<j> - <j> + <f> h e{(j> - C ) + t). (7) 
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We make the simplest assumptions about the fluctuation term £,(y,t), i.e. that it is a 
Gaussian white noise with zero mean and a correlation function given by 



where 7 denotes the nois_e s_trength. 

As was discussed in 
dimensional systems |2 



we can estimate the activation rate according to the following 
Kramers' like result for (r), the first-passage-time for the transitions between attractors, 

(n) = to exp <^ ^ } , (8) 



shown in 

As was done in 0, we assume now that the system is (adiabaticallv) subject to an external 
harmonic variation of the parameter^: 4> c (t) = 0* + 5<p c cos(ut) 0,0], and exploit the "two- 
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where AF^cj), k, L] = T[4> u {y)i k, L] — jF[0j(y), k, L] {i = 0, 1). The pre-factor r is usually 
determined by the curvature of k, L] at its extreme and typically is, in one hand, 
several orders of magnitude smaller than the average time (r), while on the other -around 
the bistable point- does not change significatively when varying the system's parameters. 
Hence, in order to simplify the analysis, we assume here that To is constant and scale it out 
of our results. The behavior of (r) as a function of the different parameters (k, 4> c /4>h) was 

ically) 

state approximation" Q as in 0,0,0). Such an approximation basically consist in reducing 
the whole dynamics on the bistable potential landscape to a one where the transitions occurs 
between two states: the ones associated to the bottom of each well, while the dynamics is 
contained onlv in the transition rates. For all details on the general two-state approximation 
we refer to [8J. 

Up to first-order in the amplitude 5<f) c (assumed to be small in order to have a sub- 
threshold periodic input) the transition rates W{ adopt the form 

.^-^{.^M} (9) 

where 

AJP% k, L, t) = AJ*[0, k, L) + 6<f> c ( dAJ7 y » k > L l ) cosM). (10) 

V d<p c / <t>c=4>* c 

This yields for the transition probabilities 

1 / 6(f) c \ 

Wi ~ -[jii^oti— cos{ut)y (11) 



where 



and 



Hi ~ exp 



ft; 



Ar[(p,k,L} 

7 



1,2). Using Eq. ©, it is clear that 



can be obtained analytically. 



These results allows us to calculate the autocorrelation function, the power spectrum 
density and finally the SNR, that we indicate by R. The details of the calculation were 



shown in |£| and will not be repeated here. For R, and up to the relevant (second) order in 



the signal amplitude 5(p c , we obtain 



R 



8] 



7T (ft 2 [l\ + ftl /i 2 )' 



4/i ly U 2 /Ul + /i 2 

Due to the form of ftj, we can reduce the previous expression to 

71 ^112 



(12) 



47 2 /ii + /i 2 



where 



L/2 



L/2 



(j>hO((j) s t(y) - C ) cfy 



[2<te(£)f 



(13) 



(14) 



We have now all the elements required to analyze the problem of SSSR. 

Figure 4 shows the typical behavior of SR, but now -in the horizontal axis- the noise 
intensity is replaced by the the system length L, for fixed values of k, 7 (the noise intensity) 
and the ratio 4> c /4>h (that in our scaled system is a single parameter). Such a response 
is the expected one for a system exhibiting SSSR. Within the context of NEP, it results 
clear that, in this kind of systems, the phenomenon arises due to the breaking of the NEP's 
potential symmetry. This means that, as shown in Fig. 2, due to the variation of L, both 
attractors can exchange their relative stability. For a value L = L*, both stable structures, 
the nonhomogeneous (f>i(y) and the trivial 0o(?/); have the same value for the NEP. When 
L < L* , 4>i{y) becomes less stable than (po(y), making the transitions from (f>i(y) to (po(y) 
"easier" (the barrier is lower) than in the reverse direction, reducing the system's response. 
When L ~ 0.72, (pi(y) and (p u (y) coalesce and disappear, and the response is strictly zero 
(within the linear response implicit in the two state approximation). When L > L*, <p\{y) 
becomes more stable than <fto(y), making now the transitions from <po(y) to <fii(y) "easier" 
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FIG. 4: SNR, vs. L, for k = 3., 7 = 0.1, D = 1., and c /0 ft = 0.193. 

than in the reverse direction, again reducing the system's response. Clearly, the system's 
response has a maximum when both attractors have the same stability (L = L*), and decays 
when departing from that situation. Hence, for this system and within this framework, SSSR 
arises as a particular case of the more general discussion done in j^. 

We can analyze the same problem from an alternative point of view. That is, studying 
the scaling of the NEP in Eq. (@J) with L. Due to the similarities of the present problem 
with the one discussed in Section IIV[ we stop here the discussion, and left such a kind of 
analysis to treat that problem. 

To conclude this section as well as for completeness, we change the point of view. In Fig. 
5 we show the curves of the SNR as a function of k, while keeping fixed values of L, and 
z. When k is not too large, indicating a high degree of reflectiveness at the boundary (that 
is, a reduced exchange with the environment), we see that the SNR changes for k varying 
from low to larger values, showing a broad resonance like curve. Remember that a large 
value of k indicates that the system boundaries become absorbent. Such a broadening of 
the resonance indicates the robustness of the systems' response when varying k, a parameter 
that somehow indicates a degree of coupling with the environment. However, according to 
the previous argument -about the breaking of NEP's symmetry- from the behavior in Fig. 
3 this is again the expected result. 
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3 4 5 6 7 

k 



FIG. 5: SNR, vs. k for L = 1.2, 7 = 0.1, D = 1., and fc/fa = 0.193. 

III. THE GLOBAL COUPLING MODEL 

In this section we consider one of the models discussed in [19] from the point of view 
of the NEP approach. The model we refer to is described by a set of (globally) coupled 
nonlinear bistable oscillators 

N 



xj = Xj - x] + — J2(x k -Xj) + y/2i £j(t) + fj(t), 



k=l 

d 



X 



[/({*},*) + V2 7 (15) 



3 ' 9xj 

with fj(t) = Acos(ut), {x} = (xi,X2, ..,ijv), £j(t) are Gaussian noises with zero mean and 
(£7 = 8ji 8(t — f), and where we have defined 

N 

U({x},t) = U Q ({x}) - Acos(ujt) y^Xj 

3=1 

N / 4 2\ N N N 



1 J 



= E (t - 1 ) + 57? E .) 2 - > ~m E 

i=l v 7 i=i fe=i 3=1 

N N N N 

= ^u Q (x j ) + — ^2^2(x k -x j f-Acos(ujt)^x j . (16) 

j=l j=l k=l j=\ 

Due to the structure of Eq. (JT5|l it is clear that U ({x}), the potential function in Eqs. 
(115116)1 . is the discrete form of the NEP for this problem. For A = the stationary distri- 
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bution of the multidimensional Fokker-Planck equation associated to Eq. (|15|) results 

p**m) « ex P y ( i7) 



7 

This potential has two attractors corresponding to X\ — X2 — ••• — %N — 

±1, and a barrier 

separating them at x\ = x 2 = ... = x^ = 0. 

Now, exploiting the same scheme as before but, as both attractors have the same "energy" , 
reduced to the symmetric case, we get for the SNR 

JtMnp (-^W)) M ^ Bq ,(-j^!W), (18) 

where X is a kind of "collective coordinate" (the one evolving along the trajectory joining 
both attractors, that pass though the saddle, and that can be approximately interpreted as 
X « J2f=i x j- However, we need the evaluation at only two points: X = 0, ±1), and 

Au (X) = u (X = ±1) - u (X = 0). 

Such a SNR clearly shows similar SSSR characteristics as those described in [19:]. In this 
situation the NEP's symmetry is retained when varying N, while the wells are deepened (or 
the barrier separating them is enhanced). However, if we scaled out N, we find a constant 
"effective" potential (uo{X)), while the system shows an effective scaling of the noise with 
N. In this case we could speak of a noise scaled SSSR, in contrast to the previous case that 
could be called a NEP symmetry breaking SSSR. 

In order to deepen our understanding of this case, let us analyze a continuous model, that 
is tightly connected with the previous discrete one. Consider a field ■?/> (?/,£), that behaves 
according to the following functional equation 

^-My,t) = ^(y,t)-^(y,t) 3 + e / &(y',t) - if>{y,t)) dy' + £(y,t) + f(t), 



dt Jn 

= -^W(y.*))+e(y,*), (19) 

with f(t) = A cos(u;t), while for the noise we assume that, as before, it is white and Gaussian 
with (£,{y,t)) = 0, and the correlation 

(ayM(y\t')) = 2 1 5(y-y')6(t-t>). 
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Q indicates the integration range, and s ^ y ^ is a functional derivative. We consider a finite 
system in the interval y G [-L/2, L/2], and assume Neumann boundary conditions. The 
form of the potential U(if)(y,t),t) results 

U(iP(y,t)) = U (i>{y,t))-F(iP{y,t)) 



dyu {iP{y,t)) + £ - dy dy' {^{y 1 \t) - ?p(y,t)) 2 - F{ijj{y,t)), 
1 Jn Jn 



where 



and 



(20) 

ij(y,t) 4 <i[>{y,ty 



F(iP(y,t)) = Acos(u;t) / ip(y,t)dy. 

n 

This potential is clearly similar to the one discussed in |7], but with the local (diffusive) 
coupling being zero, and the nonlocal contribution becoming "global". For A = the 
stationary distribution of the multidimensional Fokker-Planck equation associated to Eq. 
(HH results 

Pstat{4>stat{y)) ~ exp I 1 , (21) 

with 7 the noise intensity. This potential has two attractors corresponding to the constant 
fields ipstatiy) = ±1; and a barrier separating them at ip s tat(y) = 0. 

Hence, exploiting the same scheme as in the previous section, but reduced to the sym- 
metric case as both attractors have the same "energy" , we get for the SNR 

R - exp (.mbmm) S3 I exp (_^MMV))\ . (22) 



7 / 7 V 7 
where 

Au (i) sta t(y)) = u (ip s tat(y) = ±i) - u (i) sta t(y) = o). 

This SNR clearly shows the same SSSR characteristics as those described for the discrete 
case, where the role of iV (number of elements) is now played by L (size of the system). 

Note that Eq. (JT§j) corresponds to the continuous limit of Eq. ()15j) and the result for 
the discrete case (Eq. (118)1 ) is recovered in Eq. (|22jl for the normalized noise intensity 
Jdc = 7/ Ax (with jdc the noise intensity for the discrete case), and L = NAx. 
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To conclude this section, we refer to another case discussed in |l9j, the one corresponding 
to the Ising model. Such a case has many similarities with the case of the set of coupled 
nonlinear bistable oscillators discussed above. It can be described in a similar way to the 
case above. That means we could also find an effective potential playing the role of the NEP, 
having two attractors (corresponding to all spins up or all down), a barrier corresponding to 
a mixed state, whose high depends linearly with N (the number of spins). The final result 
will be similar to the one in Eq. (fTHj) above. 



IV. MULTIPLICATIVE NOISE CASE 



A. Brief Review of the Model 

The basic model to be considered in this section is the same one studied in |15], and 
consist of the following ensemble of nonlinear coupled oscillators, described in terms of a 
continuous field 

= ^) +/<*> + 7^ «*«)■ < 23 > 

Here £(?/, t) is again a Gaussian white noise with zero mean and correlation (£(y, t)£(y', if)) = 
2j6(y — y')6(t — if), being 7 the noise intensity. D((j)) is a field-dependent diffusivity and 
he coefficient of the noise term guarantee that fluctuation-dissipation relation is fulfilled 
. The nonlinearity f{(j>) which drives the dynamics in absence of noise is monostable, 
and we adopt a density dependent diffusion coefficient to generate a noise-induced bistable 
dynamics. In particular, we use 

D <*> = TW' (24) 

and 

f(<f>) = -0 3 + b<f>, (25) 

being Do, h and b positive constants. We will consider a finite system, limited to (-L/2 < 
x < L/2), and assume Dirichlet boundary conditions 0(±L/2) = 0. 

As we are considering the Stratonovich interpretation, the stationary solution of the 
probability P s t{4>) of the stochastic field <fi{x,i) given by Eq. (J2l?j) can be written j^J in 
terms of an effective potential 

P st (<j>) ~ eX p(-K///7), (26) 
14 



with 



V. 



effl 



L/2 f | 

dy 

-L/2 



2 K D{ ^ 



U(<f>) - A \nD(<f>) 



(27) 



(where U(<p) = Jl .D (</>') /(</>') d(j)'). Here A is a renormalized parameter, related to 7 by 
A = 7/2 Ay in a square discrete lattice, where Ay is the lattice parameter [29(. The extremes 
of V e ff — stationary noise-induced structures of the effective dynamics — can be computed 
from 

^(m^)+F eff (<f>)=0 (28) 

with an effective nonlinearity 

FefM) = m + A ^p^(^) = <w- m & - 4), (29) 

where (f)\^ depend on parameters, in particular on the renormalized noise intensity A. (We 
have found one trivial homogeneous structure = and two nonhomogeneous patterns, the 
unstable (saddle) <p u and the stable one <p s (see Ref. H))- „ 

We note that in the deterministic problem we have a monostable 30| reaction term 
(A = 0). As we increase the noise intensity, due to the noise effects, we have an effective 
nonlinear term F e ff that results to be bistable (in the interval < A < Do/(2h)) and finally 
monostable for A > Dq/{2K) (reentrance effect). We also remark that <fi = is always a root 
of F e ff (see Fig. 6), and it is an extremum of V e ff[<f)] for all values of A. In what follows we 
will call this structure 0o- As a final remark, the situation here is that the same noise that 
induces the patterns and the bistability is the one inducing the transitions among them and 
the SR phenomenon. 



B. System Size Stochastic Resonance 

To analytically describe the stochastic resonance, again we resort to use the two state 
approach in the adiabatic limit As indicated before, all details about the procedure 
and the evaluation of the SNR could be found in jSj. The system is now subject to a time 
periodic subthreshold signal b = b + S(t) where S(t) = Ab sin(a;ot). Up to first-order in 
the small amplitude Ab the transition rates W{ take the form 

W\{t) = /ix — cti Ab sin(a;o t), 
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FIG. 6: Form of the nonlinearities for the deterministic case (A = 0), bistable case (A = 0.8) 
and a monostable case (A = 1.2) in the reentrance region. The vertical scale was changed in the 
deterministic case in order to clarify the figure. Note that = remains as a root in all cases. 
The parameters used are: Dq = 1, h = 1/2 and 6 = 2. 



W 2 (t) = /i 2 + «2 sin(cu t), 



(30) 



where the constants ui ■> and ai 2 are obtained from the Kramers-like formula for the tran- 

n 

sition rates [26J 

|^g] 1/2 exp[-(^(*)-^(^))/7]. (31) 

Here /?+ is the unstable eigenvalue of the deterministic flux at the relevant saddle point (<f> B ) 
and 



2tt 



Ml,2 = Wi j2 \s(t)=0, 

dW h2 . 

\S(t)=0- 



«1,2 



dt 



(32) 



We note in passing that, due to the system's sensitivity to small variations in the param- 
eters, and at variance with the case studied in Section II, here we require the evaluation of 
the pre-factor in Eq. (pTTj) . 

As before, these results allows us to calculate the autocorrelation function, the power 
spectrum and finally the SNR (indicated by R) . For R, up to the relevant (second) order in 
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the signal amplitude Ab, similarly to Eq. (|T2"j) and (|13|) we obtain 



R 



7T (a 2 /il + «i /i 2 ) 2 7T /il/i2 



4/ii/i 2 111 + Hi 



47 2 //I + yU 2 



where now 



L/2 (y) 
-L/2 J<t> 



(33) 



(34) 



gives a simultaneous measure of the spatial coupling (through -D (</»)) and the system size 
extension (through J dy). In our previous work 2] we have found that the dependence of 
the SNR as a function of A is maximum at the symmetric situation A = A c = 0.8, where 
both stable structures (</>o and <p a ) have the same stability (V^//[^o] = Veff[4>s])- 

To analyze the system size dependence of the SR, we fix the renormalized noise intensity 
(A = A c = 0.8) and the parameters D = 1, b = 2 and h = 1/2; only change the length L 
(with fixed lattice parameter Ax). 

It is worth here remarking that in what is varied is the length of the lattice, while 
the noise intensity, the coupling, etc, are kept constant. In the present case, due to the 
characteristics of the model, we have that as L is varied, the fields change inducing the 
change of the diffusive coupling, making a strong difference with the case in [li^ . 

At this point we can just analyze the dependence of the SNR with L as was done in 
Section |H] However, as was indicated near the end of that section, we use an alternative 
for of analysis, looking at the scaling of the potential with L. We consider the following 
t r ansf or mat ions 



y -> x = f 



L > 



L ' 



the effective potential (Eq. (j2*Tj) could be written as 

'D(4>) d 12 



V eff [<f>]=L f' 2 dx{-U(<P) + \ 

J -1/2 1 1 



L dx 



, fl/2 

I — XL dxlnD((f)), 

' J-l/2 



(35) 



that finally yields (using the previous definition D((f>)/L = Di(4>)) 



V eff {<P\/L = V sc {L)= f' 2 dx{ - U{4>) + \ 

J-l/2 ( 



- AlnL>i(0)} - AlnL. (36) 
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Here it becomes apparent that this scaling yields a logarithmic length contribution to the 
scaled potential. 

The stationary solution P s t(4>) (see Eq. EEJ) of the stochastic field <fi(y,t) can be written 

P si (0) ~ exp(-V sc (L)h x (L)), (37) 

with j x (L) = j/L. 

We can also consider the scaling of the spatial factor $ (see Eq. ^ results 

,1/2 f^*) , A , /2 ,, D /" 1/2 , f, , , arctan(v / ^0 s (x))l , , 
$ = L / da:/ D(4>) (p d(f) = L —- / cfa; <U s (x) v y v ;; V . (38) 

•/-1/2 i^o " [ V/i J 

Hence, we have the dependence of the NEP as well as transitions rates (Eq. (|31jl ) and finally 
of the SNR (Eq. (}3*3*j0 . on the system length (or the number of coupled elements for discrete 
systems). To illustrate this, in Fig. 7 we show V(<p e ff) as a function of L. The behavior 
shown in this figure is analogous to the one observed in Fig. 2. Hence, we can anticipate 
the existence SSSR in this system. 

We see that for small L, small size effects increase the NEP values of the nonhomogeneous 
structures, and the uniform state results to be the most stable one. It becomes metastable 
at L c , and nonuniform patterns are the globally stable attractors for larger values of L. 
The rate transitions also reflect this fact (they are decreasing functions of L). However, 
we expect that $ — which depends on the system size — , increases with L, and due to the 
interplay of the rates and the behavior of $, a SSSR can be expected. Such a behavior 
becomes apparent in Fig. 8. 

We can make the following interpretation in terms of V SC (L), the "effective" NEP (the 
scaled form of the NEP, Eq. ([36))). and the the scaling of the noise intensity. In Fig. 9 
we depict the form of V SC (L) as function of L. It is clear that, even though it is weak, the 
dependance of V SC (L) with L still shows the change in the relative stability of the attractors, 
while we have the scaling of the noise intensity with L. Hence, we can argue that there is a 
kind of "effective entanglement" between the symmetry breaking and the noise scaling. 

V. CONCLUSIONS 

The study of SR in extended or coupled systems, motivated by both, some experimental 
results and the growing technological interest, has recently attracted considerable attention 
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FIG. 7: Nonequilibrium potential V e ff[(p s t] evaluated in the stationary structures as a function of 
the system size L. Curves correspond to: (1) stable (<f> s ), (2) homogeneous (0o) and (3) unstable 
{<p u ) patterns. Note the global stabilization of the nontrivial stable pattern for high values of L. 




FIG. 8: Signal-to-noise ratio vs. L. As indicated in the text, we fixed D Q = 1, h = 0.5, 6 = 2, and 
A = X r = 0.8. 



12|. In previous papers 



we have studied the SR phenomenon 



for the transition between two different patterns, exploiting the concept of nonequilibrium 



potential 



13 



14. 



Recent works 



have shown that several systems presents intrinsic SR-like 
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FIG. 9: V SC {L)/L evaluated on the stationary structures as a function of the system size L. Curves 
correspond to: (1) stable (<^> s ), (2) homogeneous (4>o) and (3) unstable (4> u ) patterns. These 
corresponds to the same curves as in Fig. 7. 

phenomena as the number of units, or the size of the system is varied. This phenomenon, 
called system size stochastic resonance, has also been found in a set of globally coupled units 
described by a 4 theory and even shown to arise in opinion formation models j^J. 
Such SSSR phenomenon occurs in extended systems, hence it was clearly of great interest 
to have a description of this phenomenon within the NEP framework. 

Here, we have discussed in detail two of the cases analyzed in |2l( and presented a third, 
interesting one, that corresponds to the study of the system size dependence of SR in the 
same system studied in (lfij . 

A relevant aspect that arose from these studies is that there is a kind of entrainment 
between the symmetry breaking of the NEP as described in |2jJ], together with a scaling of 
noise intensity with system size as in [hoi ]. 

In the first case we focused on a simple reaction-diffusion model with a known form of the 
NEP 2, 2s| , and -as in all three cases- considering the adiabatic limit and exploiting the 
two-state approximation, we were able to clearly quantify the system size dependence of the 
SNR. We have shown that in this case, SSSR is associated to a NEP's symmetry breaking. 
For the second case we analyzed the model of globally coupled nonlinear oscillators discussed 
in Q], and have shown that it can also be described within the NEP framework, but now 
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SSSR arises through an "effective" scaling of the noise intensity with the system size. For 
the third studied case, we have obtained the exact form of the noise-induced patterns (both 
the stable and unstable ones) as well as the analytical expression of the NEP. The interplay 
of the transition rates, that are essentially decreasing functions of L, and the behavior of $, 
that increases with L, explain the existence of a maximum in the SNR for a specific length 
of the system and a fixed noise intensity. What arose here, through an alternative form of 
analysis, is that there is a kind of entrainment between the symmetry breaking of the NEP 



as described in 



2l| , together with a scaling of noise intensity with the system size as in 



The results found in this work clearly show that the "nonequilibrium potential" (even 
if not known in detail, see for instance j^]) offers a very useful framework to analyze a 
wide spectrum of characteristics associated to SR in spatially extended or coupled systems. 



Within t 

systems 



lis framework the phenomenon of SSSR looks, as other aspects of SR in extended 
as a natural consequence of a breaking of the symmetry of the NEP or to an 



effective scaling of the noise intensity as in Q|, or could be interpreted as an entrainment 
between the two aspects. 

In addition, in the first studied case, we have seen a new form of resonant behavior 
through the variation of the coupling with the surroundings. In such a case the system's 
response to an external signal becomes more robust, that is less sensitive to the precise value 
of the albedo parameter. This fact opens new possibilities for analyzing and interpreting 
the behavior of some biological systems j^. 

As a final comment, the main difference between the first and third cases when compared 
with the second one, is that in the two former cases we have local interactions with albedo 
b.c. (covering the range from Neumann to Dirichlet b.c), while the latter has a non local 
coupling together with boundary conditions that could be assumed as Neumann. From our 
results, it is possible to argue that the effective scaling of the noise that arises in the second 
case comes from the non local interaction, and not from the b.c. To make this aspect more 
obvious, we plan to study, within the present framework, the competence between local and 
non-local spatial couplings , that arise in some multi-component models. Also, the 
consideration of more general systems with several components will allow us to analyze the 
system size dependence of SR between patterns in general activator-inhibitor-like systems. 
All these aspects will be the subject of further work. 
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